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Abstract 

Data from HERA, LEP and the Tevatron, as well as from low energy experi- 
ments are used to constrain the scale of possible electron-quark contact interactions. 
Different models are considered, including the most general one, in which all new 
couplings can vary independently. Limits on couplings and mass scales are ex- 
tracted and upper limits on possible effects to be observed in future HERA, LEP 
and Tevatron running are estimated. Total hadronic cross-section at LEP and e~p 
scattering cross-section at HERA are strongly constrained by existing data, whereas 
large cross-section deviations are still possible for Drell-Yan lepton pair production 
at the Tevatron. 



1 Introduction 



Search for "new physics" has always been one of the most exciting subjects in the field of 
particle physics. The results presented in 1997 by the HI and ZEUS [0| experiments 
at HERA electrified the physics community. Both experiments reported an excess of 
events in positron-proton Neutral Current Deep Inelastic Scattering (NC DIS) at very 
high momentum transfer scales Q^, as compared with the predictions of the Standard 
Model. Unfortunately, in spite of the significant increase in the integrated data luminosity, 
these results have not been confirmed nor contradicted]^, The effect can be just due 
to a statistical fluctuation, but can also be a first sign of some " new physics" . 

In 1998 HERA experiments started againQ to collect electron-proton data aiming for 
integrated luminosity comparable with that of the positron-proton data. The first results 
are expected soon. The aim of the presented analysis is to review experimental and 
theoretical constraints on possible signals of "new physics" at HERA and extract limits on 
a new effects to be seen in the new HERA e~p data. Limits corresponding to other present 
and future high-energy experiments are also considered. The contact interaction models, 
used as the general framework for this analysis are described in section |^. In section 
^ the relevant data from HERA, LEP, the Tevatron and other experiments are briefly 
described. Methods used to compare data with contact interaction model predictions are 
discussed in section |[ The results of analysis within different contact interaction models, 
including extracted limits on the mass scale of new interactions, are presented in section 
|5|. Predictions for the future discovery potential at HERA, as well as at LEP and the 
Tevatron are discussed in section |[ 

The analysis presented here is based on the approach suggested in [Q. Significant 
work has been done to improve the treatment of experimental data, including a proper 
interpretation of statistical and systematic errors as well as acceptance cuts and smearing. 



^Previously HERA run in clcctron-proton mode in 1992-94. 
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2 Contact Interactions 

Four-fermion contact interactions are an effective theory, which allows us to describe, in 
the most general way, possible low energy effects coming from "new physics" at much 
higher energy scales. This includes the possible existence of second generation heavy 
weak bosons, leptoquarks as well as electron and quark compositeness [|, 0. Contact 
interactions can be represented as additional terms in the Standard Model Lagrangian 0: 

Lci = Vs{eLeR){qLqR) + v'si^LeRjigRqi) + h.c. scalar (1) 

+ Vij{^i1^(^i){Qj1fi(lj) vector 

i,j=L,R 

+ riT{eL(y^''eR){qL(Tf,uqR) + h.c. tensor 

where subsequent lines describe the scalar, vector and tensor contact interaction terms 
respectively. As very strong limits have been already placed on both scalar and tensor 
terms this paper considers vector terms only. 

The influence of the vector contact interactions on the ep NC DIS cross-section can 
be described as an additional term in the tree level eg eq scattering amplitude 0: 



^ ^ ~ t ^ sin^Ow ■ cos^Qw t-Ml + "^'^ 



12 



where t = —Q is the Mandelstam variable describing the four-momentum transfer be- 
tween the electron and the quark, Cg is the electric charge of the quark in units of the 
elementary charge and the subscripts i and j label the chiralities of the initial lepton and 
quark respectively: i,j = L,R. gf and gj are electroweak couplings of the electron and 
the quark 

9l = hf - efsin^ Ow (3) 

where I^f is the third component of the SU(2) isospin for the fermion /: / = e, q. 

For processes such as e~^e~ hadrons or pp l^l~X, a corresponding formula can 
be written for the e^e~ — > qq tree level amplitude: 

= _4vr«e^ ^ . T#^-TV + Vi] (4) 

s siri^Ow ■ cos^Ow s - + isj^ ^ 
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where s is the center-of-mass energy squared of the four-fermion reaction. The sign of 
the contact interaction contribution to the s-channel amphtude (Q) is the same as for the 
i-channel amphtude (Q). However, the Standard Model amphtude changes its sign due 
to the opposite signs of s and t variables. It is therefore important to notice that the 
resulting sign of the interference terms in the cross-section for e^p scattering is different 
from that in e~^e~ or pp scattering. 

The contact interaction coupling strength rj can be related to the mass scaleQ Ai of 
new physics through the formula: 



V 




where gci is the unknown coupling strength of new interactions. As the contact interaction 
contribution always depends on the gci to M. ratio, it is convenient to consider the 
effective mass scale A defined through the formula: 

47r 

which corresponds to the choice g^j = Att. 



2.1 General Model 

In the most general case, vector contact interactions are described by 4 independent cou- 
plings for every lepton-quark pair. With only 2 lepton (e and fi) and 5 quark flavours (i.e. 
neglecting t quark contribution), we still have 40 independent couplings. It would be very 
difficult, if not impossible, to consider the model with 40 free parameters. However, some 
of these parameters (couplings) are weakly constrained by existing experimental data. To 
reduce the number of the free model parameters, weakly constrained couplings can be 
either neglected or additionally constrained by relating them to some other couplings. 

Most of existing experimental data is sensitive predominantly to electron-up and 
electron-down quark couplings. Therefore, the first model considered in this analysis 
is the one assuming that these 8 couplings {rjffi, Tif^R-, v'rli V'rri Vll^ vTri V'rli V'rr) can 
vary independently, whereas other couplings (for s, c, 6, t quarks and/or yU, r leptons) are 
assumed to vanish. This case will be referred to as the general model. 
^Exchanged particle mass or compositeness scale. 
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The other possibihty is to impose additional relations between couplings. The common 
choice is to assume lepton universality: 



and quark family universality: 



C (6) 



Lepton universality allows us to include data on muon pair production at the Tevatron 
(see section p.2|) , whereas assuming quark family universality significantly improves the 
constraints which we can obtain from LEP2 measurements (see section p.3| ). As a result, 
experimental constraints on contact interactions can be significantly improved without 
increasing the number of free model parameters. The model assuming relations (^ and 
d^) will be referred to as the model with family universality. 



2.2 SU{2)l X U{1)y Universality 

Another commonly used assumption about lepton-quark contact interactions is that they 
satisfy the SU{2)l x f/(l)y gauge invariance of the Standard Model. Assuming that left- 
handed electrons and quarks belong to SU{2)l doublets and that the contact interaction 
Lagrangian (|lD respects the SU{2)l symmetry implies a relation between contact terms 
involving left-handed u and d quarks [§: 

^eii _ ed 

'Irl — 'Irl 

which reduces the number of free model parameters from 8 to 7. The SU{2)l x U{1)y 
also relates eeqq contact interaction couplings with those of uuqq interactions 

vll = vli (7) 
Vli = vTl 
'Ilr — 'Ilr 

vli = vIr 
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This allows us to use, in the study of eeqq contact interactions, additional data on NC 
neutrino scattering (see section p.4| ). 

Moreover, assuming the SU{2)l x U{X)y universality introduces a related contact 
interaction term in the Charged Current process eg — > vq' . The coupling constant for the 
induced Charged Current contact interaction is 

^cc^^ew ^ r^l\-,^^^^ (8) 

This relation allows us to use, in the study of Neutral Current contact interaction, also 
data from Charged Current processes (see section [O] and |3^ . 

The model assuming the SU {2)^ x U{1)y universality will be referred to as the SU(2) 
model. In order to reduce the number of models, the SU(2) models considered in this 
analysis always assume lepton and quark family universality. 



2.3 One-parameter models 

Using data from a single experiment it is mostly not possible to put significant constraints 
on contact interaction scales in the general case. Therefore it is a common practise to 
consider particular models, which assume fixed relations between the separate couplings, 
reducing the number of free parameters to one. For example, the so called vector- vector 
model assumes that all couplings are equal: 

Mass scale limits obtained in one-parameter models are, artificially, much stronger than 
in the general model. They will be considered in this analysis to allow comparison with 
other results. The relations between couplings assumed for different models are listed in 
Table It should be noticed that all one-parameter models considered assume 

eq I eq eq eq n 

Vll + Vlr-Vrl-Vrr = 
for q = u,d, to avoid strong limits coming from atomic parity violation measurements 



(see section |3.4|) . For all one-parameter models quark and lepton family universality 
is assumed. The results obtained both with and without imposing the SU{2)l x U{1)y 
universality are presented, except for the U2, U4 and U6 models, which violate it explicitly 

(Vrl^Vrl)- 



5 



Global Analysis of eeqq Contact Interactions and Future Prospects for High Energy Physics 



"I\ T 11 
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Table 1: Relations between couplings for the one-parameter models considered in this 
paper. 

3 Experimental Data 

In this section the data used to constrain contact interaction model are presented. For 
each measurement, the formula describing the possible influence of the new couplings on 
the measured quantities is given. Description of the statistical methods used to interpret 
the data will be presented in the section ^j. 



3.1 High-Q2 DIS at HERA 

Used in this analysis are the latest data on high-Q^ e~^p NC DIS from both HI and 
ZEUS IQ, corresponding to integrated data luminosities of 37 and 47p6^\ respectively. 
Older results from e~p NC DIS scattering ^ are also used, although the influence of 
these data is marginal. For models with the SU{2)l x f/(l)y universality, as mentioned 
in section 0, data on e^p CC DIS 0, are also included in the fit. 



HERA experiments quote their high-Q DIS results in terms of numbers of events 
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and/or cross- sect ionsQ measured in bins of Q^. For simplicity let us consider a single 
bin ranging from (5^j„ to Q^ax- Assume that usm events are expected from the Standard 
Model. 

The leading order doubly-differential cross-section for positron-proton NC DIS {e~^p 
e'^X) can be written as p 

{mil' + |M^5,r + (1 - y? [mA\' + mil'] } 

where x is the Bjorken variable, describing the fraction of proton momentum carried by 
a quark (antiquark), q{x) and q{x) are the quark and antiquark momentum distribution 
functions in the proton and M^j' are the scattering amplitudes of equation which can 
include contributions from contact interactions described by a set of couplings ff. 

The cross-section (including the contribution from contact interactions), integrated 
over the x and range of an experimental bin is 

Qmax 1 j2 LO ( ~\ 



Q 



2 

min 



where i/max is the upper limit on reconstructed Bjorken variable y imposed in the analy- 
sis0. The number of events expected from the Standard Model with contact interaction 
contributions can now be calculated as: 

n{v) = nsM-[^^] (10) 

where cx^^j is the Standard Model cross-section calculated with formula (P) (setting r/ = 
0). Leading-order expectations of the contact interaction models are used to rescale 
the Standard Model prediction nsu coming from detailed experiment simulation. This 
accounts not only for different experimental effects, but also for higher order QCD and 
electroweak corrections. Validity of this approach is discussed in section ^ 

^If not given, the number of events can be estimated from the cross-section value assuming that the 

statistical error quoted corresponds to the Poisson error on the number of measured events, ctjv — ^fN ■ 
'^For the NC DIS analysis HI uses Umax = 0.9, whereas ZEUS uses Umax — 0.95. For CC DIS analysis 

both experiments use ymax — 0.9. 



7 



Global Analysis of eeqq Contact Interactions and Future Prospects for High Energy Physics 



3.2 Drell-Yan lepton pair production at the Tevatron 



Used in this analysis are data on Drell-Yan lepton pair production from the CDF ||T2 



and DO |]T3[ experiments. Both experiments present numbers of measured high-mass 
electron pairs {pp — > e^e~ X). CDF also presents results on muon pair production 
{pp — > fJ'^fi'^ X), which are used in the case of models with family universality (see 
section H). 

The leading order cross-section for lepton pair production in pp collisions is 

where Mu is the invariant mass of lepton pair, Y is the rapidity of the lepton pair center- 
of-mass frame, Xi and X2 are the fractions of proton and antiproton momenta carried by 
the annihilating quarks. The scattering amplitudes M^j and the parton density functions 
are calculated for scale 

S = X1X2S 

where s is the total proton-antiproton center of mass energy squared. 

The cross-section corresponding to the Mu range from Mmm to M^ax is calculated as 

.-(,7) ^ / I "VMY}-^ (11) 

where Ymax is the upper limit on the rapidity of the produced lepton pair: 

Y = In^ 

max iJ^-"- , r ; 

and Aii{Y) is the acceptance function, resulting from the integration over the lepton pair 
production angle in the center of mass system, with angular detector coverage taken into 
account. The cross-section calculated with equation ([TT|) is used to calculate the number 
of events expected from the Standard Model with contact interaction contributions using 
formula (|T0|). 

3.3 Measurements from LEP 

Many measurements at LEP are sensitive to different kinds of "new physics". The eeqq 
contact interactions can be directly tested in the measurement of the total cross-section 



8 



A.F.Zarnecki 



3 EXPERIMENTAL DATA 



for e"'"e~ — >■ qq. Using flavour tagging tecliniques, additional constraints can be obtained 
from the measurement of the heavy quark decay fractions Rb and Rc, and of the forward- 
backward asymmetries A'p^ of qq events. 

The leading order formula for the total quark pair production cross-section e+e~ qq, 
at the total electron-positron center of mass energy squared s, is 

where M^j' are the scattering amphtudes described by equation (H), including contribu- 
tions from contact interaction couphngs rf. For comparison with measured experimental 
values, the leading order contact interaction cross-sections are rescaled using the expected 
Standard Model cross-section a^^\s) quoted by experiments: 

where cr^*^(s, 0) is the leading-order Standard Model cross-section (f/ = 0), calculated with 
equation (p!^). This takes into account possible experimental effects and higher order QCD 
and electroweak corrections (for discussion see section^). All four LEP experiments have 



recently presented data on ahad for center-of-mass energies up to 189 GeV |Tj, |T^, [l6|, |17 



The sensitivity of the total hadronic cross-section to the contact interaction coupling 
strength ff is limited by the fact that the interference terms in the quark-pair production 
cross-sections have opposite signs for up-type and down-type quarks. In the total cross- 
section, summed over all quark flavours]^, these terms tend to compensate each other. 
However, if this is the case, the fraction of events produced with the given quark-pair 
flavour turns out to be very sensitive to the contact interaction couplings. 

Using different flavour tagging techniques, cross-sections for bb and cc pair production 
and the corresponding fractions Rb and Rc can be measured. Although the limited tagging 
efficiency and purity significantly affects the measurement, useful constraints on contact 
interaction couplings can be extracted. Used in this analysis are results on Rb coming 
from Aleph[0, |T^, Delphi ||15| and Opal fg^j as well as Delphi results on Rc [p!5[. 



'Production of the t quark is taken into account only for y/s > 350 GeV. 
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The contact interaction contribution to the scattering amphtude affects also the ob- 
served forward-backward asymmetry of qq events. In the leading order the forward- 
backward asymmetry can be calculated as 



mi 




Mil 


2 _ 






2 


\Mll\^ + \Mll 


2 + 




2 + 




2 



where the factor | corresponds to the integration over the full angular rangeQ. 

Constraints upon the forward-backward asymmetries A'p^ are obtained using a jet 
charge technique. After clustering all the events into two jets, the jet charge Qjet of each 
jet can be determined from the momentum weighted sum over all charged tracks in the jet. 
The sign of Qjet coincides with the charge of produced quark in about 70% of events. The 
forward-backward asymmetry for the selected sample of events (e.g. 6-tagged events) can 
be extracted in two ways. The method used by Aleph is based on the measurement of the 
mean charge difference between the forward and backward jets {Qfb) = {Qfet) " (Qfet)- 
Delphi and Opal extract A'^p^ from the angular distribution of jets with well defined 
sign. In both cases, the measured asymmetry depends on the parton-level asymmetries 
A'p^ and on the quark content of the selected sample. As the up-type and down-type 
quarks have charges of opposite signs, the measured asymmetry is very sensitive to the 
relative contribution of different quark flavours. Even if we measure the asymmetry for 
the flavour-tagged sample, the selected sample of events is always contaminated by other 
quark flavours (e.g. a 6-tagged sample always contains a fraction of cc events) and the 
measured value depends strongly on the quark production fractions (e.g. Rb and Rc). 
This is the reason why the measurement of the forward-backward asymmetry is very 
sensitive to the contact interaction couplings. Used in this analysis are the measurements 
of forward-backward asymmetry for the 6-tagged events |^ 0, c-tagged events ||T5[| 



and anti-tagged events |14, Q 



3.4 Data from low energy experiments 

The low energy data are included in the present analysis in the manner which follows 



closely the approach presented in g, ^ . Therefore only basic assumptions are listed here 
^For results which are based on the sample of events selected with | cos 6*1 < 0.9, this factor is reduced 
to 0.70866... 
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and technical details are omitted. 

In case of the general contact interaction model the following constraints from low 
energy experiments are considered: 

• Atomic Parity Violation (APV) 

The Standard Model predicts parity non-conservation in atoms caused (in lowest 
order) by the Z° exchange between electrons and quarks in the nucleus. Experimen- 
tal results on parity violation in atoms are given in terms of the weak charge Qw 
of the nuclei. A very precise determination of Qw for Cesium atoms was recently 



reported |22|. The experimental result differs from the Standard Model prediction 
H, a by: 

AQ^' = Q^"""' - Q^ = 0.71 ±0.84 



Corresponding results have also been obtained for thallium |25|, |24 

AQ^ = 1.9 ±3.6. 

These measurements are used to place limits on contact interaction contributions 
to Qw- 

97 -I- 

+ ^J2G ^^^^ ^ ^'^^ ~ ^'^^ ~ ^rr) 

electron-nucleus scattering 

The limits on possible contact interaction contributions to electron-nucleus scatter- 
ing at low energies can be extracted from the polarisation asymmetry measurement 

^ _ daR — d(TL 
daR + dai 

where da^^R^ denotes the differential cross-section of left- (right-) handed electron 
scattering. Polarisation asymmetry directly measures the parity violation resulting 
from the interference between the weak {Z° exchange) and the electro-magnetic (7 
exchange) scattering amplitudes. For isoscalar targets, taking into account valence 
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quark contributions only, the polarisation asymmetry for elastic electron scattering 
is 



A. 



20na, 



2 {gl + 9r) - {gi + gi 



where is the four-momentum transfer and gf are quark electroweak couplings, 
as introduced in equation (Q). Contact interactions modify the effective quark elec- 
troweak coupling 



gi 



gl 



eg 
VLi 



2V2G1 



(14) 



The constraints used in this analysis come from the SLAC eD experiment [26|, the 
Bates eC experiment and the Mainz experiment on eBe scattering In case 
of models with family universality also data from the /x^C experiment at CERN|]29[| 
are included]]]. 

In case of the SU(2) models additional constraints come from: 

• neutrino-nucleus scattering 

Constraints on the couplings of quarks to the Z° and/ or additional uugq contact 
interactions (related to eeqq CI, as described in section can also be derived 
from the precise measurement of the ratio of Neutral Current to Charged Current 
neutrino-nucleon scattering cross sections 



a 



NC 



a 



uN 

CC 



However, when using constraints on gf resulting from measurement of R", one also 
has to take into account that possible contact interaction contribution affects not 
only the Neutral Current but also the Charged Current scattering cross-section (see 
section |2.2D . Therefore the quark electroweak coupling extracted from R'^ measure- 
ments should be expressed asPI 



gi 



gl 



uq 



1 - 



nCC 



4V2Gp 

^The constraints from the /U^C experiment result from the comparison of n^N and n]^N cross sections. 
^This correction seems to be missing in M, |21|. 
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It is important to notice that 77^^ entering formula (|T4D has been replaced here by rj^j. 
This is because the effective g1 and gf couplings measured in neutrino scattering are 
sensitive to "flavour crossed" contact interaction couplings r^f^^ and yyf^ respectively, 
which results from relations (J^). Experimental constraints on i?^ come mainly from 
muon-neutrino experiments. Assuming lepton and quark family universality, the 
following measurements of gf from z/^A^ scattering are used: the results compiled 



by Fogi and Haidt |30| and the recent constraints from CCFR|31] and NuTeVp2| 



lepton-hadron universality of weak Charged Currents 

Charged Current contact interactions which are induced by SU{2)l x U{1)y uni- 
versality (see equation (^) would also affect the measurement of Vud element of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix, leading to the effective violation of 
unitarity |3^, |3^. The current experimental constraint is ||2^ 



\Vud? + \Vus? + \Vub? = 0.9969 ± 0.0022 



whereas the expected contribution from the contact interaction is 

Vud 



VuT-(i 



aV2G 



electron-muon universality 

In the similar way Charged Current contact interactions would also lead to effective 



violation of e-fi universality in charged pion decay The current experimental 
value of R = T{Tr~ — > ez/)/r(7r~ fiu) is ||3^ 



R 



R. 



SM ■ 



= 0.9966 ±0.030 

RsM 

whereas the expected contribution from the contact interaction is 

^cc \ 2 

It is interesting to notice, that data in Charged Current sector may point to a slight 
violations of the unitarity of the CKM matrix and of the e-fi universality. Both measure- 
ments are consistent with the presence of CC contact interactions with a mass scale of 
the order of 10 TeV. The combined significance of these two results is about 1.8a, but it 
has a considerable influence on global analysis results for the SU(2) model. 
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4 Analysis method 

3.4 Data from low energy experiments 

The aim of this study is to find the allowed ranges for contact interaction couplings 
within the different models considered. To do so, the probability function in the coupling 
space, 

-Piv) ~ I[P^iv), (15) 

i 

is calculated. In ([T5|) , the product runs over all experimental data i and ff represents the 
set of free parameters for a given model (one or many). This section describes how the 
probability function is defined and which corrections are included to take into account 
experimental conditions. 

4.1 Statistical errors 

All experimental data used in this analysis can be divided into two classes. 

1. For experiments in which a result can be presented as a single number with an error 
which is considered to reflect a Gaussian probability distribution, the constraints 
on the contact interaction couplings can be usually expressed using the equation 

F{if} = AA ± aA 

where AA is the difference between the measured value and the Standard Model pre- 
diction, and F{ff) is the expected contact interaction contribution to the measured 
value of A. The resulting probability function can be written as 



1 (F(r/) - AAf 

2 ^ 



- expi--^ '-'"-. 1 (16) 



reflecting the definition of the Gaussian error a a- This approach is used for all low 
energy data as well as for the LEP hadronic cross-section measurements. 

2. On the other hand, when the experimentally measured quantity is the number of 
events of a particular kind (e.g. HERA high-Q^ events or Drell-Yan lepton pairs at 
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the Tevatron), and especially when this number is small, the probability is better 
described by the Poisson distribution 

p.(^) _ n(r;)^ ■ exp{-n{f])) ^^^^ 

where N and n(f/) are the measured and expected number of events in a given 
experiment, respectively, and n(r/) takes into account a possible contact interaction 
contribution. This approach has been used for HERA and the Tevatron data. 

4.2 Systematic errors 

For low energy data the total measurement error can be used in formula (|TBp taking into 
account both statistical and systematic errors. For collider data, formula (|T6|) or ( [T7| ) 
is used to take into account the statistical error of the measurement only. As for the 
systematic errors, it is assumed that within a given data set (e.g. e~^p NC DIS data from 
ZEUS ) they are correlated to 100%. This seems to be a much better approximation of the 
experimental conditions than assuming that systematic errors are uncorrelated[|. In fact 
most of the contributing systematic uncertainties at HERA are highly correlated between 
different bins, as for example energy scale uncertainty or the luminosity measurement. 

For each data set, a common systematic shift parameter 6 has been introduced to 
describe the possible variation of event numbers expected at HERA or the Tevatron, or 
cross-sections predicted at LEP, due to systematic error: 

nsM = nsM + 
or a''' = asM + 5 ■ 

where fish[ (o'sm) is the nominal expectation from the Standard Model and o"^^** (c'ct^'*) is 
the total systematic uncertainty attributed to this number. Parameters S can been treated 
as additional free parameters when maximising the overall model probability V{rf). When 
doing so, normal probability distributions for parameters 5 are included in the definition 
(|l5|) of the probability function]^. 



^Unfortunately the experiments do not publish the correlation matrix for their systematic errors so 

these are the only possible choices. 

^''This corresponds to the assumption that systematic errors are described by the Gaussian probability 

distribution 
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4.3 Migration corrections 

Equation (|r^) introduced in section ^ takes into account experimental conditions at HERA 
and the Tevatron. The number of events expected with contact interaction contribution 
is calculated by rescaling the Standard Model prediction usm coming from the detailed 
experiment simulation. However, this is only an approximation based on the assumption 
that the acceptance for contact interaction events is the same as for standard NC DIS 
events. Although the detection efficiency for given or {Mii,Y) is always the same 

(as we have the same final state), the distribution of events in the kinematic plane in the 
presence of the contact interactions can differ significantly. This can affect the measure- 
ment due to the finite or Mu resolution. To take this effect into account a dedicated 
migration correction is introduced. The DIS cross-section in the bin from to 
Qmax is calculated by the following extension of eq. (0): 



S 1 

^msiv) = JdQ^-S{Q^;Ql,^,Ql,,,aQ2) J dx 



d^ajfj) 
dxdQ"^ 



where (Jq2 is the resolution, as quoted by experiments, assumed to be constant within 
the bin. The Drell-Yan cross-section is calculated by the similar extension of eq. (|TTp: 

^dy{v) = J dMu ■ S{Mu; M^,„, M^,,, a,,) J dY MY) ■ 

-Y 

^ max 

where gm is the Mu resolution. The mass resolution has been estimated from the quoted 
calorimeter energy resolution (for electrons) or tracking momentum resolution (for muon 
pairs). The acceptance function used in both formula 



S{x\ a, b,a)= J dy 



2710- 



l(y-a)2\ / liy-b)^' 



exp — — — exp 



2 ; ^ V 2 a2 



describes the probability that the true value x measured with resolution a will be re- 
constructed between a and b. The migration corrections are important for the muon-pair 
production results from the Tevatron and for the CC DIS results from HERA. For electron- 
pair production or for NC DIS results, when the corresponding mass and resolutions 
are much better, the effects of the migration corrections are very small. 



16 



A.F.Zarnecki 



4 ANALYSIS METHOD 



The influence of tlie systematic errors and tlie introduced smearing on tlie model 
probability function V{ff) has been studied for the ZEUS e^p NC DIS data [Q. The 
results, in terms of the log-likelihood function, -logP, for four chosen one-parameter 
models, are shown in Figure |l|. The applied corrections (mainly the systematic error 
correction) can have sizable influence on the model probability distribution. Taking into 
account statistical errors only leads to much narrower probability distribution and gives 
much stronger constraints. The most prominent effect is observed for the VV model. A 
narrow probability maximum (minimum of -logV function) observed when only statistical 
errors are included, becomes wider with a "shoulder" on one side when the systematic 
errors are taken into account. 

The results from this analysis have been compared with the ZEUS results based on 
full detector simulation The comparison for the same four one-parameter models is 
presented in Figure 0. For some models very good agreement is observed between this 
analysis and ZEUS results, as can be seen for AA and XI models. However, for models 
such as VV or U2, the constraints given by ZEUS are stronger (probability distribution 
narrower) than the constraints resulting from this analysis. This is due to the fact that 
the ZEUS analysis takes into account the two-dimensional event distribution in the {x, y) 
plane, whereas this analysis uses the one- dimensional distribution only. 

4.4 Radiative corrections 

For high-energy data from HERA, LEP and the Tevatron, Standard Model predictions 
given by experiments are used to rescale leading-order expectations of the contact interac- 
tion models (see formula ( p!OD and (|13])). This accounts not only for different experimental 
effects, but also for higher order QCD and electroweak corrections, including radiative cor- 
rections. This approach is reasonable as long as the difference between the corrections for 
the Standard Model and for the model including contact interactions is negligible. It is 
natural to assume that this difference should be much smaller than the correction itself. 

The contribution of radiative corrections to high-Q^ DIS at HERA is of the order of 
10%. For high-mass Drell-Yan lepton pair production at the Tevatron it is only about 
6%. Therefore, the possible variation of the radiative corrections for both HERA and 
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Figure 1: Log-likelihood function -Xo'gVij]) for ZEUS e+p NC DIS data, for four chosen 
one-parameter models, as indicated on the plot. The functions are calculated with statis- 
tical errors only (dashed line) and with migration and systematic error corrections (solid 
hne) . 
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Figure 2: Log-likelihood function -log 7^(77) for ZEUS e+p NC DIS data, for four chosen 
one-parameter models, as indicated on the plot. The results from this analysis (solid line) 
are compared with the ZEUS results obtained with full detector simulation (dashed line). 



19 



Global Analysis of eeqq Contact Interactions and Future Prospects for High Energy Physics 



Tevatron data have been neglected. The only data where radiative corrections could be 
significant is the hadronic cross-section measurement at LEP. 

Most of the events observed at LEP2 are radiative events. This is due to the "radiative 
escape" to the Z° peak. Radiation probability is significantly enhanced as the e^e~ 
annihilation cross-section at ^/s = Mz is several orders of magnitude higher than at 
nominal i/i. The leading order cross-section (^) is corrected for radiation effects using 



the formula p6| 



s s 



where integration runs over the center-of-mass energy squared s' of the produced quark 
pair, and s^j„ is the minimum value of s' required by the event selection cuts[5 G(z) 
is the "radiator function" encapsulating the results of QED virtual and real corrections. 
Used in this analysis is the approximate formula (based on |3^, |3^ ) 

Giz) = f^,.p{l-zy~' + il-fr)-6il-z) 

where /3 = 2 i log-^ — 1 ) 

TT \ mi J 

The parameter fr is chosen to reproduce the cross-section ratio for radiative and non- 
radiative eventsQ. 

It turned out that the effect of radiative corrections on the probability function V{ff) 
is very small. The resulting limits on contact interaction mass parameters decrease by at 
most 3%. 



4.5 Probability functions 

The probability function V{ff) summarises our current experimental knowledge about 

possible eeqq contact interactions. It will be used to set limits on contact interaction 

mass scale parameters and to extract predictions concerning possible future discoveries. 

It is therefore very important to understand the precise meaning of Vlfj). 

"Data used in this analysis correspond to ^Js'/s > 0.9 (ALEPH) or y^s'/s > 0.85 (DELPHI, L3 and 
OPAL). This choice significantly reduces possible influence of radiative corrections. 
i^Events with O.K y^TJl <0.85 and y^VJl >0.85 
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V{'q) is not a probability distribution of f/. A probability distribution should describe 
the probability of finding a given value of variable. Our situation is different. Function 
Vijf) describes the probability that our data come from the model described by the set 
of couplings f] (see section [4.1|) . It is our data set, which is a variable, and r/ is a set of 
model parameters: they are unknown, but they are fixed. This simple observation has 
very important implications for this analysis, not only for the limit setting procedure (see 
next subsection) but also for calculation of model predictions. 

To set limits on possible deviations from the Standard Model predictions (eg. for 
NC DIS cross-section at very high-Q^ at HERA or for hadronic cross-section at next 
e^e~ collider), we have to consider the probability function P(r), where the cross-section 
deviation r is defined as 

r = ^ = P(r/) 

<^SM 

If V{f]) is taken as probability distribution, then the probability distribution for r should 
be calculated as 



P(r) = j d^T] V{f]) 6{r - R{f])) 



where integration is performed over N-dimensional coupling space. This however leads 
to completely false results, as is demonstrated in appendix Instead of calculating the 
probability distribution for r (which is not well defined) , we should rather try to find out 
what is the probability that our data come from the model predicting deviation r. This 
leads to the formula: 

P(r) = {vm 

R(ff)=r 

where averaging is necessary, if we want to reduce number of parameters of the probability 
function (for multi-parameter models). The commonly used assumption in that case, 
is that 7] has fiat underlying (prior) distribution^. The formula for P(r) can be then 
expressed as 

/ d^ffV\v) 6{r~R{ff)) 

P{r) = (19) 

d^vV{i]) 6{r-R{r])) 



^■^This corresponds to the assumption, that we would have no preferences for any vahie of 77, if there is 
no experimental data. 
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The formula applies for any variable which can be used as a parameter of the probability 
function. In this analysis it will also be used to calculate probability functions and to 
set limits on mass scale parameters corresponding to single couplings in multi-parameter 
models. 

As V{r]) is not the probability distribution it does not satisfy any normalisation con- 
dition. Instead it is convenient to rescale the probability function in such a way that its 
global maximum has the value of 1: 

maxV{ff) = 1. (20) 

4.6 Extracting limits 

After imposing condition (pO]), the lower and upper limits on the value of the model 
parameter r are defined as minimum (r~) and maximum {r~^ ) values satisfying relation 

P(r-) = 0.05 
and P(r+) = 0.05 . 

For any model described by the parameter r < r~ or r > r+, the probability that our 
data results from this model is less than 5% of the maximum probability. This is taken 
as the definition of the 95% confidence level (CL) limits. 

For one-parameter contact interaction models this approach is slightly modified. As 
models with negative and positive values of t] are usually considered as independent 
scenarios (differing by the signs of the interference terms in the cross-section), the upper 
and lower limits on rj are calculated using restricted rj range: 

p(rj-) = 0.05-maxP(r7) 

and P(7]+) = 0.05 ■ max P(77) . (21) 

»7>0 

For one-parameter contact interaction models, or for probability functions related to 
single couplings in multi-parameter models, the limits on coupling values t]~ and t]'^ can 
be translated into the limits on contact interaction mass scales 



A- 
A- 



4tt 
A-n 
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Mass limits commonly used in literature are based on t] limits defined in a slightly 
different way. In this paper they will be denoted as rj and rj^^. Their definition follows 
from the equations: 



J dr] P{ri) = 0.95 ■ J drj P{r]) 

rj-- -oo 

oo 

J dr] P{r]) = 0.95 • J dr] P{r]) (22) 



This approach is based on the assumption that r] has a flat underlying (prior) distribution. 
In such a case P{r]) can be treated as the probability distribution for rj. The mass scale 
limits corresponding to rj and r]'^^ will be denoted as A and A^"*". Although the def- 
inition resulting from equation (plf ) is considered to be more appropriate for this analysis 
than definition (^2]) , the results for both definitions are presented to allow comparison 
with other results. 

As definitions ( pT]) and (p2D correspond to the different interpretation of the probability 
function, they are not expected to give similar results. In fact, the allowed range for 
parameter rj, calculated with equation (^) is usually about 25% wider than the one 
calculated with equation (^)Q. As a result, corresponding mass scale limits A^ and A+ 
are usually 10 to 15% smaller than A and A++. 



5 Results 

For one-parameter models the analysis has been performed both without and with the 
additional SU{2)l x U{1)y universality assumption (see section ||). In the latter case 
data coming from neutrino-nucleus scattering experiments and from different Charged 
Current processes ( refer section ^ have been also used to constraint contact interaction 
couplings. In the following, the models assuming SU{2)l x U{1)y symmetry are referred 
to as SU(2) models. One-parameter models without SU{2)l x U{1)y symmetry will be 
referred to as simple models, to avoid possible confusion. For all one-parameter models 

quark and lepton family universality is assumed. 

^''For Gaussian shape of the probabihty function, rj^ and 77+ correspond to ±2.45(j Umits, whereas rj 
and r]++ to ±1.96ct. 
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Mo del 


V 


Mass scale limits 


[ieVJ 






A- 


A+ 


A— 


A++ 


VV 


-0.015 ±0.049 


9.8 


10.7 


11.0 


11.8 


AA 


0.007 ±0.048 


10.5 


10.1 


11.7 


11.3 


VA 


0.049 ±0.143 


6.6 


6.2 


7.3 


6.9 


XI 


0.014 ±0.073 


8 7 


8 1 


9 6 


9.2 


X2 


-0.011 ±0.075 


8 2 


8 4 


9 2 


9.4 


X3 


-0.003 ±0.051 


9 9 


10.2 


11.1 


11.4 


X4 


-0.113 ±0.138 


0. ( 


K 9 

o.z 


0.41: 


5.3 


X5 


-0.079 ±0.132 


5.9 


6.4 


6.6 


7.0 


X6 


-0.013 ±0.147 


6.2 


5.8 


7.0 


6.4 


Ul 


-0.059 ±0.104 


6.4 


7.7 


7.3 


8.4 


U2 


-0.065 ±0.082 


6.9 


9.1 


7.8 


9.9 


U3 


-0.044 ±0.053 


8.5 


11.7 


9.6 


12.7 


U4 


-0.136 ±0.166 


5.1 


5.5 


5.7 


5.8 


U5 


-0.093 ±0.092 


6.4 


8.8 


7.2 


9.5 


U6 


0.115 ±0.128 


7.0 


5.6 


7.4 


6.3 



Table 2: Coupling values and 95% CL mass scale limits resulting from fits of one-parameter 
models without SU{2)l x U{1)y universality. The errors attributed to rj values corre- 
spond to the decrease in the model probability V{ri) by the factor of ^/e. See text for 
explanation of the symbols. 

Using the overall model probability P(r/), as defined by equation ([15|), the "best" 
values of contact interaction couplings (i.e. corresponding to the maximum probability) 
were found using the MINUIT package The results for one-parameter simple and 

SU(2) models are presented in Tables ^ and ^ respectively. The errors attributed to rj 
values correspond to the decrease in the model probability V{ri) by the factor of y/e. In 
case of asymmetric errors the arithmetic mean is given. 

For all simple one-parameter models considered couplings are found to be consistent 
with the Standard Model within la. The same is true for most SU(2) models. However, 
for the SU(2) models Ul and U3 the "best" coupling values are more than 2cr from the 
Standard Model. These are the only two models which allow for rj^'^ ^ (i.e. r/f^ 7^ 
r^l'f^). The observed deviation from the Standard Model predictions is directly related to 
r/'"'" bounds coming from the unitarity of the CKM matrix and the e-/i universality, as 
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Mo del 


V 


Mass scale limits 


ieVJ 






A- 


A+ 


A— 


A++ 


VV 


-0.024 ±0.047 


9.6 


11.4 


10.8 


12.5 


AA 


-0.010 ±0.047 


9.9 


11.1 


11.1 


12.3 


VA 


-0.078 ±0.108 


6.3 


8.0 


7.1 


8.7 


XI 


-0.025 ±0.067 


O.L 


Q K 


Q 9 


10.5 


X2 


-0.041 ±0.069 


7.8 


9.6 


8.8 


10.5 


X3 


-0.019 ±0.049 


9.5 


11.1 


10.7 


12.2 


X4 


-0.066 ±0.144 


6.0 


5.4 


6.7 


5.8 


X5 


-0.040 ±0.131 


6.2 


6.4 


7.0 


7.0 


X6 


-0.013 ±0.147 


6.2 


5.8 


7.0 


6.4 


Ul 


-0.100 ±0.042 


7.9 


17.0 


8.6 


17.8 


U3 


-0.083 ±0.036 


8.6 


18.2 


9.4 


19.1 


U5 


-0.050 ±0.082 


7.1 


8.8 


8.0 


9.6 



Table 3: Coupling values and mass scale limits resulting from fits of one-parameter models 
with SU{2)l X U{1)y universality. The errors attributed to rj values correspond to the 
decrease in the model probability V{ri) by the factor of ^/e. See text for explanation of 
the symbols. 

described in section ^. However, it has to be noticed that other data also do support 
this effect: the discrepancy observed for the combined data is more significant than for 
the Charged Current sector only. Although the effect is interesting, the data are still in 
acceptable agreement with the Standard Model. The probability that our data result from 
the Standard Model equals 5.7% and 7.0% for the Ul and U3 SU(2) models respectively. 
Assuming that there is no direct evidence for eeqq contact interactions, the limits on 
single couplings can be calculated. The results on mass scale limits A~, A+, A and 
A^"*" obtained from fitting one-parameter models are summarised in Tables ^ and ^. For 
simple models mass limits range from 5.1 TeV (A~ for the U4 model) to 11.7 TeV (A"^ 
for the U3 model). Similar limits are obtained for most of the SU(2) models. Only for 
the Ul and U3 SU(2) models much higher A+ limits of 17.0 and 18.2 TeV are obtained. 
The probability functions V{r]) for four selected SU(2) models are shown in Figure ^. 

Contribution of different data sets to the mass scale limits presented in Tables ^ and 
^ can be estimated using the probability function. Mass scale limits A~ and A+, derived 
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Figure 3: Probability functions V{rj) for chosen one-parameter models with SU{2)l x 
C/(l)y universality, as indicated on the plot. 



26 



A.F.Zarnecki 



5 RESULTS 



from coupling limits ?7~ and 77"*", correspond to the decrease of the global probability V{ri) 
to 0.05 of its maximum value (see equation (pT])). This decrease can be presented as a 
product of contributions from all data samples. Table ^ presents the relative probability 
changes, calculated separately for different data sets, corresponding to mass scale limits 
A~ and A"*", for different one-parameter SU(2) models. The product of numbers in every 
row is equal to the factor 0.05 defining the 95% CL. Numbers close to 1.0 demonstrate 
that given data set has negligible influence on the considered mass scale limit. The smaller 
the number, the more sensitive are the data to a given CI model. Numbers greater than 
1.0 indicate that the model with mass scale A^ or A+ gives better description of given 
data set than the "best" coupling value resulting from the combined fit. 

The results presented in Table ^ show that in most models the strongest constraints 
on contact interaction couplings come from LEP data on forward-backward asymmetries 
A'^p^ and on quark production ratios Rq. However, for particular models a significant 
contribution can result from LEP hadronic cross-section measurements, neutrino-nucleus 
scattering data, HERA NC DIS data or from data on Charged Current interactions. 

The probability functions for single couplings obtained for the general model are pre- 
sented in Figure All couplings are consistent with the Standard Model predictions. 
Results for single couplings obtained for the general model, the model with family univer- 
sality and the SU(2) model with family universality are summarised in Table ^. It has to 
be stressed that all limits for single couplings are derived without any assumptions con- 
cerning the remaining couplings, which corresponds to the definition (0) of a probability 
function. For this reason most calculated limits are weaker than in case of one-parameter 
models. The mass limits obtained for the general model range from 2.1 TeV (A^i+) to 5.1 
TeV (A^*^"). For the SU(2) model with family universality, the corresponding numbers 
are 3.5 and 7.8 TeV for A^'^^" and A^^"^, respectively. 

Single couplings can either increase or decrease cross-section for a given process, as 
compared with the Standard Model expectations. It is therefore also possible that the 
influence of the two different couplings compensate each other. Because of that, the limits 
on mass scales obtained for single couplings does not exclude contact interactions with 
smaller mass scales. To obtain the most general limit, the eigenvectors of the correla- 
tion matrix (obtained from MINUIT from the functional form of V{ff) in the vicinity of 
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Mass Relative change in model probability 



Model 


scale 


HERA 


Tevatron 




LEP 




Low energy NC 


CC 




limit 


e p iNU 


Urell- Yan 


'^had 






l^N 


vN 


data 


VV 


A- 

A -1- 

A^^ 


0.697 
1.107 


0.887 
0.632 


1.344 

0.751 


1.665 
0.205 


0.028 
0.841 


1.000 
1.000 


1.300 
0.553 




AA 


A- 
A+ 


1.240 
0.774 


1.021 

0.835 


0.536 
1.451 


1.867 
0.208 


0.020 
0.722 


0.984 
1.013 


1.961 
0.350 




VA 


A- 

A -1- 

A+ 


1.264 
1.008 


0.710 
0.928 


1.281 
1.370 


1.092 
0.811 


0.033 
0.504 


0.916 
1.059 


1.316 
0.090 




XI 


A- 
A+ 


1.214 
0.830 


1.013 
0.885 


0.694 
1.395 


1.746 
0.344 


0.017 
0.634 


0.962 
1.031 


2.060 
0.217 




X2 


A- 

A 4- 

A+ 


0.842 
1.079 


0.931 
0.748 


1.323 
0.878 


1.659 
0.327 


0.017 
0.660 


0.973 
1.021 


1.773 
0.320 




X3 


A- 

A -U 

A+ 


1.016 
0.977 


1.005 
0.710 


0.993 
1.137 


1.734 
0.192 


0.017 
0.738 


0.992 
1.007 


1.720 
0.444 






A- 
A+ 


0.331 
0.814 


0.762 
0.666 


0.318 
0.218 


1.021 
0.779 


1.414 
0.690 


1.018 
0.969 


0.424 
0.813 




X5 


A- 
A+ 


0.617 
1.142 


0.650 
0.516 


1.178 
0.884 


1.749 
0.064 


0.133 
1.411 


1.043 
0.949 


0.435 
1.114 




X6 


A- 
A+ 


0.771 
1.367 


0.786 
0.731 


1.127 
0.789 


0.077 
1.976 


0.984 
0.031 


0.965 
1.025 






Ul 


A- 
A+ 


1.141 
0.933 


0.963 
0.954 


0.895 
0.753 


1.207 
0.903 


0.361 
1.174 


0.925 
1.028 


0.443 
0.323 


0.285 
0.212 


US 


A- 
A+ 


1.056 
0.963 


0.855 
0.885 


0.254 
0.502 


1.289 
0.804 


0.243 
1.212 


0.987 
1.005 


1.157 
0.449 


0.611 
0.266 


U5 


A- 

A+ 


0.724 
1.004 


0.849 
0.671 


0.447 
0.109 


1.535 
0.570 


0.561 
1.267 


1.084 
0.910 


0.195 
1.036 





Table 4: Relative changes in model probabilities calculated for separate data sets (as 
indicated in the table) corresponding to the decrease in the global probability at mass 
scale limit (A~ or A"^) to 0.05 of its maximum value (for negative or positive couplings 
respectively). Considered are one-parameter models with family and SU{2)l x U{1)y 
universality. 
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Figure 4: Probability functions P{r)) for single contact interaction couplings (as indicated 
on the plot) obtained within the general contact interaction model. 
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Mass scale limits [TeV] 



Model with SU(2) model with 

Coupling General model family universality family universality 



A- A+ A— A++ A- A+ A— A++ A' A+ A" A++ 

rili 3.1 3.6 3.4 4.0 4.3 5.4 4.8 6.0 7.5 6.1 8.1 6.8 

rjfji 2.4 2.2 2.6 2.5 2.8 3.2 3.0 3.6 3.6 3.7 4.0 4.1 

r/lf^ 2.6 2.1 2.8 2.4 3.9 2.8 4.4 3.1 4.5 3.9 5.0 4.4 

r]'j(ji 2.8 2.8 3.0 3.1 3.3 3.5 3.7 3.9 3.5 5.1 3.8 5.7 

rjll 5.1 3.9 5.6 4.3 5.3 4.1 5.9 4.5 6.5 7.8 7.3 8.5 

r/i'lj 4.1 3.3 4.3 3.7 3.7 3.5 4.1 3.9 4.9 4.7 5.4 5.3 

rjl^L 3.1 3.5 3.5 3.9 3.1 3.5 3.4 3.9 4.5 3.9 5.0 4.4 

r^g); 3.7 3.7 4.1 4.1 3.8 4.2 4.3 4.6 4.8 4.4 5.4 4.8 

Table 5: 95% CL mass scale limits for single couplings, obtained within different models, 
as indicated in the table. See text for mass scale limits definition. 

the maximum probability) are considered. In case of the general model the two least 
constrained linear coupling combinations are 

r/i = -0.267^21 + 0.84r/2^^ + 0.15r7|i + 0.33r/|'^^ 
-0.06Cl + OAIvTr - O.lOr/- + 0.27r/- 
and ri2 = +0.20?72'1 + 0.17r/2^^ + 0.81//^^ - 0.48r/|f^ 

+0.10Cl - 0.14Cii + - 

This is in agreement with the observation that the least constrained single couplings are 
Tjf^j^ and rjl^i, which can also be seen from Table ^. The probability functions for r]i and 
rj2 are shown in Figure ^ The mass scale limit corresponding to rji is[^ 

Ai = 2.1 TeV 

This limit should be considered to be the most general one, as it is valid for any combi- 
nation of couplings. This means that any contact interaction with a mass scale below 2.1 
TeV is excluded on 95% CL. On the other hand it also shows that the existing data do 
not exclude mass scales of the order of 3 TeV. The limits on the mass scale associated 
with rji are summarised in Table ^ 



15 



As the sign of rji is arbitrary, only one value is given. It is calculated as min(A^, Aj^ 
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-0.04 -0.02 0.02 0.04 -0.04 -0.02 0.02 0.04 



Ti^pv /47t 11^^^/471 [TeV^] 

Figure 5: Probability functions, calculated within the general contact interaction model, 
for the two least constrained coupling combinations rji and 772 (upper plots), the atomic 
parity violating combination rjAPv (lower left) and the CC contact interaction coupling 
^cc induced in the SU(2) model (lower right plot). Note different horizontal scales be- 
tween upper and lower plots. 
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Mass scale limits [TeV] 

Model with SU(2) model with 

Coupling General model family universality family universality 

A- A+ A— A++ A- A+ A"- A++ A" A+ A" A++ 

r^i 2.1 2.3 2.7 3.0 3.1 3.5 

T]APv 9.8 6.0 10.4 6.6 9.7 6.1 10.3 6.7 10.9 7.5 11.7 8.4 
Vli-Vll'^-^ 3.5 3.1 3.9 3.4 4.3 3.8 4.8 14.4 7.2 15.1 7.9 

Table 6: 95% CL mass-scale limits corresponding to the least constrained coupling combi- 
nation rji, atomic parity violating coupling combination tjapv and rjf^^^ — Vll combination 
corresponding to Charged Current contact interaction coupling r]'^'^ of SU(2) model. As 
the sign of t]i is arbitrary, only one value is given. 

Shown in the same table are mass limits corresponding to the atomic parity violating 
combination of couplings, 



VAPV ^ vli + vfR-V%i-V%'R 



+ vII + vTr-Vrl-Vrr 



rjAPv is close to the most strongly constrained coupling combination (the eigenvector 
with the highest eigenvalue). Mass scale limits up to about 11 TeV are obtained. The 
probability function for rjAPv in the case of the general model is included in Figure ^ 

Also shown in Figure ^ is the probability function for the Charged Current contact 
interaction coupling ri^'~^ induced in the SU(2) model. The discrepancy between the data 
and the Standard Model has decreased slightly, as compared with the Ul and U3 SU(2) 
models. The most probable value of r/*"*" is about 2a from the Standard Model value, 
which corresponds to the probability of about 10%. This discrepancy is observed for the 
SU(2) model only. When SU{2)l x f/(l)y universality is not assumed (i.e. in case of 
the general model and the model with family universality) the corresponding coupling 
combination is no longer related to the Charged Current sector and is in good agreement 
with the Standard Model. The corresponding mass scale limits are included in Table 0. 
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6 Predictions 

All presented results are in good agreement with the Standard Model. Nevertheless, an 
interesting question is whether "new physics" in terms of contact interactions can be 
expected to show up in high-energy experiments in the near future. 

The cross-sections corresponding to the "best fit" of the general model (the set of 
coupling values resulting in the best description of all data, i.e. corresponding to the 
maximum probability) are compared in Figure ^ with the HERA, LEP and Tevatron 
data. In the case of LEP data, the best fit of the general model agrees very well with 
the Standard Model. The Contact Interaction contribution to the measured cross-section 
does not exceed 3% for ^/s up to 200 GeV. On the other hand, the same model predicts 
for both HERA and the Tevatron an increase in the cross-section by almost a factor of 
2 at the highest Q'^/Mu. In order to verify the significance of these predictions it is 
unavoidable to consider the statistical uncertainty of these predictions. 

Employing the Monte Carlo techniques, the probability function for the contact inter- 
action couplings, Vlfj), is translated into the probability function for relevant cross-section 
deviations, as described in section Considered in this analysis are possible deviations 
from the Standard Model predictions for high-Q^ e~p and e~^p scattering at HERA[^ (see 
section |3.1j ), for the total quark pair production cross-section at LEP (or Next Linear 
Collider, NLC; see section ^.3| ) and for the Drell-Yan lepton pair production at the Teva- 
tron (see section p^) . The probability functions calculated for these processes at selected 
energy scales are presented in Figure ^ 

The results for HERA, in terms of the 95% confidence limit bands on the ratio of 
predicted and the Standard Model cross-sections as a function of Q^, are shown in Figures 
P and 1^ for the general model and the SU(2) model with family universality, respectively. 
For the e~^p NC DIS the uncertainty of these predictions is very big, although the nominal 
predictions of both models are above the Standard Model. The Standard Model prediction 
is well within the 95% confidence level band. For the general model, the increase in the 
e+p NC DIS cross-section at HERA by up to about 80% at of 30,000 GeV^ would still 
be consistent with current experimental data. For the SU(2) model the corresponding 

^^For the proton beam energy of 920 GeV and the electron/positron beam energy of 27.5 GeV. 
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NC DIS at HERA 
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Figure 6: Cross-section deviations from the Standard Model resulting from the general 
contact interaction model fit (thick solid line) compared with HERA, LEP and the Teva- 
tron data. 
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1 2 3 1 2 3 

G(e*p^e*X)/asM a(e"p^e"X)/asM 




a(e*e"^qq)/asjfl o(pp^n")/OsM 

Figure 7: Probability functions for possible deviations from the Standard Model predic- 
tions for: e~^p and e~p NC DIS cross-section at HERA, at = 30,000 GeV^ (upper 
plots), e+e" total hadronic cross-section at ^/s — 400 GeV (lower left plot) and Drell-Yan 
lepton pair production cross-section at the Tevatron, at Mu — 500 GeV (lower right plot). 
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Figure 8: The 95% CL limit band on the ratio of predicted to the Standard Model cross- 
section for e'^p and e~p NC DIS scattering at HERA (upper plots) and the 68% and 95% 
CL contours for the possible deviations for scattering of right- and left-handed electrons 
and positrons at =30,000 GeV^ (lower plots). The limits are calculated using the 
general contact interaction model. 
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Figure 9: The 95% CL limit band on the ratio of predicted to the Standard Model cross- 
section for e^p and e~p NC DIS scattering at HERA (upper plots) and the 68% and 95% 
CL contours for the possible deviations for scattering of right- and left-handed electrons 
and positrons at =30,000 GeV^ (lower plots). The hmits are calculated using the 
SU(2) contact interaction model with family universality. 
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limit is 63%. It turns out that the best statistical sensitivity (in single measurement) 
to possible contact interaction effects is obtained when considering the number of events 
measured for > 15,000 GeV^. The allowed increase in the integrated e~^p NC DIS 
cross-section is about 40% for the general model and about 30% for the SU(2) model. In 
order to reach the level of statistical precision, which would allow them to confirm possible 
discrepancy of this siz^^ HERA experiments would have to collect e^p luminosities of 
the order of 100-200 pb~^ (depending on the model). This will be possible after the HERA 
upgrade planned for year 2000. 

Constraints on the possible deviations from the Standard Model predictions are much 
stronger in case of e~p NC DIS. This is because the Standard Model cross-section itself is 
higher, and also because different contact interaction coupling combinations contribute. 
It is interesting to notice that the possible cross-section increase for e~^p NC DIS, which is 
suggested by global fit results, corresponds to decrease in the NC DIS cross-section for 
e~p. For the general model, deviations larger than about 20% are excluded for > 15,000 
GeV^, whereas for the SU(2) model with family universality the limit goes down to about 
7%. When compared with the predicted statistical precision of the future HERA data, 
this indicates that it will be very hard to detect contact interactions in the future HERA 
e~p running. For the general model the required luminosity is of the order of 400 pb^^. 

However, the HERA " discovery window" can be visibly enlarged if we consider scatter- 
ing of polarised electrons and/or positrons. The 68% and 95% CL contours for the allowed 
deviations for scattering of right- and left-handed electrons or positrons are included in 
Figures || (for general model) and | (for SU(2) model), at Q"^ =30,000 GeV^. In both 
cases, the cross-section deviations for ej^p and e~^p scattering are less constrained than 
in case of e~^p and ej^p, respectively. For the general model possible deviations for both 
left- and right-handed projectiles are significantly higher than in the unpolarised case. 
However, for the SU(2) model, constraints significantly weaker than in the unpolarised 
case are obtained only for e^p and e^p scattering. In both models deviations of up to 
about 50% are still allowed for e~j^p scattering at Q'^ >15,000 GeV^, assuming 60% polar- 

"'^^We require that the allowed increase in the cross-section for > 15,000 GeV^ (at 95% CL) should 
correspond to at least three times the statistical error on the number of events. 5% systematic uncertainty 
on the expected number of events is assumed. 



38 



A.F.Zarnecki 



6 PREDICTIONS 



isation of the positron beam. To observe effects of this size it would be enough to collect 
luminosity of the order of 70-80 pb~^. For e]^p scattering maximum allowed deviations 
are 28% and 19%, for the general and SU(2) models respectively. It means that with 60% 
longitudinal polarisation it would be possible to observe significant deviations from 
Standard Model predictions already for luminosities of the order of 120 pb^^ (for the 
general model). Unfortunately, polarisation can result in significantly higher systematic 
uncertainties of the Standard Model predictions, which was not considered here. 

Since the only visible inconsistency between data and the Standard Model is observed 
in the Charged Current sector (for models assuming SU{2)l x f/(l)y universality), the 
interesting question is whether any effect can be observed in high CC DIS at HERA. 
It turns out that the possible effect is far beyond the HERA sensitivity. The "best" 
rj'-^^ value (resulting from the SU(2) model fit) corresponds to a decrease in the CC DIS 
cross-section at HERA not greater than 2% within the accessible range, and a decrease 
exceeding 5% is excluded at 95% CL. At the same confidence level, any increase in the 
cross-section by a similar amount is excluded. 

Model predictions for both the total hadronic cross-section at electron-positron collider 
(LEP or NLC) and Drell-Yan lepton pair production cross-section at the Tevatron are 



shown in Figures |T0| and |TT], for the general contact interaction model and the SU(2) 
model, respectively. For e~e'^ qq at y/s above about 300 GeV upper cross-section limits 
obtained from both contact interaction models increase rapidly. Cross-section deviations 
up to a factor of 3 are allowed for ^/s ~500 GeV. Unfortunately, this energy range will 
be accessible only in the Next Linear Collider experiment (s). LEP will not go beyond 
^/s ~200 GeV: at this energy the possible deviations from the Standard Model are only 
about 8%, which makes possible discovery very difficult. 

However, significant deviations from Standard Model predictions are still possible for 
heavy quark production ratios Rc and Rb, and for the forward-backward asymmetries 
Ap^ and A'^p^. The 68% and 95% CL contours for the values of the forward-backward 
asymmetry versus the quark production fraction, allowed within the SU(2) model for c 
and b quark production at y^=200 GeV are included in Figure |1T[ Allowed ranges for Rc, 
Rb, Ap^ and A^p^ at y^=200 GeV, for different contact interaction models considered, are 
summarised in Table ^. For the general model, in which contact interactions are limited 
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Figure 10: Left: the 95% CL limit band on the ratio <Jsm+ci/o'sm, where asM is the 
Standard Model total hadronic cross-section at LEP/NLC and ctsm+ci is the cross-section 
calculated in the general contact interaction model; Right: the same ratio for the Drell- 
Yan lepton pair production at the Tevatron. 





SM 


Allowed range on 95% CL. 




Value 


General 


Model with 


SU(2) model 




(LO) 


model 


family univ. 


w. fam. univ. 


Rb 


0.159 


0.147 -0.161 


0.137-0.180 


0.139 -0.179 


Rc 


0.262 


0.242 -0.266 


0.230 -0.294 


0.232 -0.291 


^FB 


0.601 




0.345 -0.750 


0.431 -0.732 


^FB 


0.668 




0.469 -0.750 


0.551 -0.738 



Table 7: Leading order Standard Model prediction and the allowed range (on 95% CL) 
for the heavy quark production ratios and forward-backward asymmetries, for e~^e~ an- 
nihilation at ^/s=200 GeV, in different contact interaction models. 
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Figure 11: Upper left: the 95% CL limit band on the ratio o'sm+ci/(^sm, where asM is the 
Standard Model total hadronic cross-section at LEP/NLC and ctsm+ci is the cross-section 
calculated in the SU(2) contact interaction model; Upper right: the same ratio for the 
Drell-Yan lepton pair production at the Tevatron; Below: the 68% and 95% CL contours 
for the allowed values of the forward-backward asymmetries and the quark production 
fractions for c and b quark production at LEP, at y/s=200 GeV. 
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to the first quark generation only, variations of Rb and Rc are still possible, due to the 
possible changes in uu and dd production cross-sections. However, parton level forward- 
backward asymmetries Ap^ and A^p^ do not depend on contact interaction couplings in 
this model. Therefore limits for Ap^ and A^p^ are not reported for the general model. 
Heavy quark observables considered here are least constrained for the model with family 
universality. Large effects are still possible in this model for both production fractions and 
asymmetries. Deviations up to about 13% are possible for Ri, and Rc- Least constrained by 
the existing experimental data is the forward-backward asymmetry for the bb production 
A''p^, where deviations from the Standard Model prediction by up to 40% are still possible. 
Note that for the general model Rb and Rc are 100% correlated, whereas for models with 
the family universality they are 100% anti-correlated. 

It seems that the best place to study contact interactions in the nearest future is the 
Tevatron, which should run again after being upgraded in the year 2000. If there is any 
"new physics" corresponding to the contact interaction model it is very likely to show 
up in Drell-Yan lepton pair production for masses above 200-300 GeV. Moreover, upper 
limits on possible deviations from the Standard Model predictions are much higher than in 
case of HERA and LEP/NLC. For Mu =500 GeV, which should be easily accessible with 
increased luminosity, cross-section deviations up to a factor of 5 are still not excluded. 

Upper limits on the cross-section deviations from the Standard Model predictions, 
derived on 95% confidence level in different contact interaction models are summarised in 
Table |. 

When considering possible future discoveries at high-energy experiments, it is also 
interesting to study the relation between effects observed at different experiments. The 
68% and 95% CL contours for the sizes of the allowed deviation from the Standard Model 
predictions, for different measurement combinations, are shown in Figures [1^ and 0, for 
the general contact interaction model and for the SU(2) model with family universality, 
respectively. In both cases, clear correlation is observed between the Drell-Yan cross- 
section deviation at the Tevatron and the hadronic e~^e~ cross-section at LEP/NLC. 
Possible cross-section increase at the Tevatron has to be accompanied by the increase 
in the hadronic cross-section at LEP/NLC. Similar correlation is observed between the 
hadronic e~^e~ cross-section at LEP/NLC and e~^p NC DIS cross-section at HERA for 
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Limits on Aa/asM [%] 



Reaction 


Energy 


General 


Model with 


SU(2) model 




scale 


model 


family univ. 


w. family univ. 


e+p NC DIS 


(5^=10000 GeV^ 


11 


10 


9 


^/s =318 GeV 


Q2=20000 GeV^ 


36 


30 


28 




Q2=30000 GcV^ 


81 


65 


63 




g2=50000 GcV^ 


220 


180 


170 


e-p NC DIS 


Q2=ioOOO GeV2 


8 


4 


3 


=318 GeV 


g2=20000 GeV^ 


18 


8 


7 




g2=30000 GeV^ 


28 


13 


11 




Q2=50000 GeV^ 


49 


26 


21 


e~e'^ — > qq 


v/i=175 GeV 


5 


5 


6 




v^=200 GcV 


8 


8 


8 




Vi=225 GeV 


14 


13 


11 




^/s=250 GeV 


26 


24 


16 




v/i=300 GeV 


65 


61 


35 




v^=400 GeV 


185 


185 


110 


pp i+i-x 


M,z=200 GeV 


17 


12 


12 


Vi = 1800 GeV 


M,i=300 GeV 


64 


55 


38 




M;,=400 GeV 


190 


185 


95 




Mh=500 GeV 


440 


450 


210 



Tabic 8: Upper limits (on 95% CL) on cross-section deviations from the Standard Model 
predictions in different contact interaction models. Considered are e^p and e~p scattering 
at HERA, total hadronic cross-section at LEP/NLC and Drell-Yan lepton pair production 
at the Tevatron, as indicated in the table. 
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Figure 12: The 68% and 95% CL contours for the possible deviation from the Standard 
Model predictions, for different combinations of measurements. Considered are: e~^p and 
e~p NC DIS cross-section at HERA, at = 30,000 GeV^, total e~^e^ qq cross-section 
at ^/s = 400 GeV and Drell-Yan lepton pair production cross-section at the Tevatron, 
at Mil = 500 GeV, as indicated on the plot. The limits are calculated using the general 
contact interaction model. 



44 



A.F.Zarnecki 



6 PREDICTIONS 



m 
D 



O 68% CL 
n 95% CL 



0) 



1.5 



3, 



0.5 - 




0.8 1 1.2 1.4 

o(e"p^e"X)/o, 



SM 



0.6 0.8 1 1.2 1.4 

o(e"p^e"X)/o, 



SM 



0.8 1 1.2 1.4 

o(e"p^e'X)/o, 



12 3 4 

o(e*e"^qq)/Os„ 




1 1.5 2 

o(e>^e*X)/o, 



SM 




SM 



1 1.5 2 

o(e>^e*X)/o, 



SM 



Figure 13: The 68% and 95% CL contours for the possible deviation from the Standard 
Model predictions, for different combinations of measurements. Considered are: e^p and 
e~p NC DIS cross-section at HERA, at = 30,000 GeV^, total e+e~ qq cross-section 
at ^/s = 400 GeV and Drell-Yan lepton pair production cross-section at the Tevatron, 
at Ma = 500 GeV, as indicated on the plot. The limits are calculated using the SU(2) 
contact interaction model with family universality. 
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SU(2) model. Another interesting observation is that the possible decrease in the e'p 
NC DIS cross-section at HERA should be related to the increase in both Tevatron and 
LEP/NLC cross-section. In other cases correlations between different measurements are 
weak. This shows that contact interaction searches at LEP, the Tevatron and HERA 
are, to large extent, independent. Data from all types of experiments are necessary to 
constraint contact interaction model in general case. 
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7 Summary 

Data from HERA, LEP, the Tevatron and low energy experiments were used to constrain 
electron-quark contact interactions. The contact interaction mass scale limits obtained 
for different one-parameter models range from 5.1 to about 18 TeV. Using the most gen- 
eral approach, in which all couplings to are allowed to vary independently, any contact 
interactions with mass scale below 2.1 TeV are excluded at 95% CL. This limit can be 
raised to 3.1 TeV by assuming SU{2)l x U{1)y and quark/lepton family universahty. 
There is a shght hint on possible "new physics" in the Charged Current sector (related 
to Neutral Current contact interactions by SU{2)i x f/(l)y universality), where the dis- 
crepancy between data and the Standard Model is at the 2(7 level. The mass scale of new 
Charged Current interactions suggested by the data is of the order of 10 TeV. However, 
this effect - if real - would have negligible impact on predictions for future collider results. 

The limits on possible effects to be observed in future HERA, LEP and Tevatron 
running are estimated. Possible deviations from the Standard-Model predictions for total 
hadronic cross-section at LEP and e~p scattering cross-section at HERA, are already 
strongly hmited by existing data. However, improved experimental sensitivity to new 
interactions should result from the measurement of heavy quark production ratios and 
asymmetries at LEP, as well as from polarised electron scattering at HERA. Sizable effects 
are still not excluded for e^p NC DIS at HERA and the required statistical precision of the 
data should be accessible after HERA upgrade. The best " discovery potential" seems to 
come from future Tevatron running, where significant deviations from the Standard Model 
predictions are still allowed. For Drell-Yan lepton pair production cross-section deviations 
at M//=500 GeV up to a factor of 5 are still not excluded. However, all experiments should 
continue to analyse their data in terms of possible new electron-quark interactions, as 
constraints resulting from different experiments are, to large extent, complementary. 
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A Interpretation of the probability function 

In this appendix a simple "toy model" is used to demonstrate that the probability function, 
as introduced in section ^, should not be treated as the probability distribution for fj. 

Let us consider a model with N independent couplings. Assume that all data con- 
sidered in the analysis are in perfect agreement with the Standard Model and that the 
resulting probability function is 



where fj"^ = rjf + ... + rjfj and the distribution width a is taken to be the same for all 
couphngs. The Standard Model gives the best description of the data, corresponding to 
the maximum value of V{'q). 

Consider the cross-section deviation from the Standard Model prediction, which is of 
the form 

r[ff) = ^ = (l+r/2) 

If V{ff) is taken as a probability distribution, then the probability distribution for r should 
be calculated from equation ([18|) . After integrating over the coupling space we obtain: 

The shape of P{r) corresponds to that of the distribution for N degrees of freedom. 
For N < 2, P{r) has a maximum for r = 1, i.e. for the Standard Model expectation. 
However, for models with > 3 parameters, the maximum of P(r) is shifted towards 
r > 1 and the probability of the Standard Model solution P{r = 1) = 0. This results is 
incompatible with our initial assumption that all data are in perfect agreement with the 
Standard Model. 

The above calculation, based on the formula (0), is not correct because it assumes 
that V{ff) is the probability distribution for rf. We can treat Vlfj) as the probability 
distribution only if we assume that ff has a flat prior distribution. This assumption 
justifies limit setting procedure described in section |4]^ (formula (0)). However, it does 
not justify variable transformation from -ff to r (resulting from equation (pISD), as the prior 
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distribution for the new variable does not need to be flat. Instead, one should try to 
define P{r) in the same way as V{f]), i.e. as the probability that our data come from the 
model predicting deviation r. This approach results in formula (|19|). For our toy model 
the probability of observing deviations from the Standard Model predictions is 

where the normalisation condition (|20|) has been imposed. The result does not depend on 
the number of free model parameters and the most probable model is the one predicting 
no deviation from the Standard Model (taking into account that r > 1). 
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